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基于 言 退 卷 积 的 AIA 图 像 增强 方法 


摘 要 : 太阳 图 像 中 存在 着 各 种 尺度 ， 亮 度 ， 结 构 不 同 的 物理 活动 现象 ， 由 于 太阳 日 网 高 动态 
的 活动 和 传感器 设备 等 因素 ,造成 太阳 图 像 成 像 质 量 不 佳 。 根据 太阳 动力 学 天 文 台 (Solar Dynamic 
Observatory, SDO) 的 大 气 成 像 仪 《Atmospheric Imaging Assenbly, AIA) 拍摄 不 同 波段 数据 结构 的 
动态 范围 大 、 骂 声 大 、 结 构 相 对 模糊 等 特点 ， 本 文 提出 一 种 有 效 基于 盲 退 卷 积 的 SDO/AIA 图 像 增 
强 方法 。 首 先 对 SDO/AIA 图 像 进行 去 噪 和 降低 动态 范围 的 处 理 ， 基 于 图 像 功率 谱 的 分 布 假设 ， 从 
原 图 中 估计 出 点 扩散 函数 (point spread function，PSE) 的 功率 谱 ， 然 后 使 用 相位 提取 算法 恢复 出 PSF 
的 相位 ， 再 退 卷 积 出 较 高 质量 的 目标 图 像 ; 最 后 通过 轮廓 切片 分 析 、 功 率 谱 分 析 以 及 PSF 分 析 对 
增强 结果 进行 定量 和 定性 分 析 。 实 验 结 果 表 明 ， 相 比 现 有 图 像 增强 方法 ， 该 方法 在 有 效 增 强 太 阳 日 
园 图 像 细 节 结 构 的 同时 ， 复 原 出 原 图 因 模 糊 无 法 辨识 的 结构 。 
关键 词 : 太阳 图 像 ; SDO/AIA ;功率 谱 ; 盲 退 卷 积 ; 相位 提取 
中 图 分 类 号 : P182.6+2 ”文献 标识 码 : A 


构 边 缘 模糊 


免 图 像 增强 。 


噪声 较 大 ， 使 得 可 视 化 困难 


斯 滤波 后 , 


Histogram Equalization, BPDHE) 这 些 方法 能 够 有 效 提高 图 像 整 体 对 比 度 。 但 日 


为 了 研究 高 度 动态 的 太阳 日 免 结构 ， 人 们 一 直 在 提高 时 间 、 空 间 分 辨 率 上 ， 不 断 探 索 。 
日 蝎 结 构 的 密度 、 位 置 和 形状 短 至 几 秒 钟 就 会 发 生变 化 ， 较 大 的 动态 范围 ， 造 成 望远镜 拍摄 的 日 狗 结 


于 太阳 


mm 


。 目 前 已 经 研究 出 儿 种 较为 有 效 的 增强 图 像 算法 用 于 太阳 
基于 直方 图 均衡 的 方法 ， 可 以 利用 不 同 分 布 直 方 图 对 不 同 图 像 规 定 化 处 理 门 ， 也 可 以 高 
再 独立 均衡 各 个 分 区 ， 如 保持 亮度 的 动态 直方 图 均衡 化 (Brightness Preserving Dynamic 


于 直方 图 的 领域 形状 


通常 是 正方 形 ， 会 导致 严重 的 伪 影 ， 从 而 影响 日 锡 环 状 结构 ， 其 次 会 放大 图 像 中 低 对 比 度 噪声 ， 进 一 


步 模 糊 图 像 细 节 ; 亮度 不 同 的 


区 域 之 间 也 会 失去 有 效 结构 。 基 于 曝光 融合 框架 的 对 比 度 增 强 算法 站 


(Exposure Fusion Framework, EFF) 可 以 解决 直方 图 均衡 出 现 的 对 比 度 过 高 或 者 过 低 的 问题 。 基 于 多 
种 经 典 图 像 增强 算法 融合 的 方法 外 (Multi-deviation Fusion method, MF) 融合 多 种 经 典 图 像 增强 算 法 


进行 取 长 补 


短 ， 使 图 像 增强 时 能 够 同时 录 


顾 细节 增强 、 对 比 度 提升 和 主观 感知 效果 等 方面 。 这 些 方法 


大 部 分 预 处 理 都 采用 对 数 变换 ， 能 够 使 得 计算 复杂 度 降低 。 但 由 于 对 数 变换 抑制 了 亮 区 梯度 幅度 的 变 


化 ， 这 可 能 会 导 


致 某 些 区 域 的 精细 结构 细节 丢失 。 为 解决 这 个 问题 ， 基 于 Retinex 的 方法 ， 可 以 压缩 动态 范 
围 压 缩 的 同时 ， 兼 顾 边缘 增强 和 保持 颜色 稳定 ， 其 中 (Robust Retinex Model , rR) 算法 握 可 提高 对 噪声 的 


鲁 棒 性 ， 使 用 一 种 基于 拉 格 朗 日 乘 子 的 交 蔡 方向 最 小 化 算法 代替 对 数 变 换 ， 有 效 解决 模型 优化 过 程 问题 ， 
能 有 效 保 持 细节 结构 和 边缘 稳定 。 上 述 方法 都 仅 为 太阳 日 冕 图 像 提 供 了 一 种 可 视 化 方法 ， 对 于 其 中 因为 


日 晃动 态 变 化 、 望 远 镜 抖动 、CCD 欠 采 样 以 及 杂 散 光 等 因素 造成 精确 结构 的 模糊 问题 , 无 法 有 效 增 强 。 
随 着 数字 图 像 处 理 技术 的 发 展 ， 盲 退 卷 积 算法 对 于 要 求 恢复 结构 准确 性 和 真实 性 方面 提供 了 科学 的 理 


论 依据 。Paolo G 等 人 中 根据 望远镜 设计 理论 估计 出 理论 上 满足 SDO/AIA 不 同 波段 的 PSF 用 于 退 卷 积 ， 


= 


的 PSF 并 不 能 很 好 的 退 卷 积 出 清晰 图 像 。 在 使 用 空间 望远镜 SDO/AIA 观测 太阳 


晶 造 成 图 像 模糊 的 原因 除 传感器 设备 外 ， 还 有 太阳 活动 位 移 、 形 状 变化 等 造成 的 模糊 。 所 以 理论 估计 


HEME, ERARE, CCD 欠 采 样 以 及 杂 散 光 等 因素 造成 成 像 模 糊 过 程 ， 可 建 模 为 : 


其 中 ，* 表 示 卷 积 运算 ， 好 表示 获得 的 模糊 图 


B(x) =1(x)*k(x)+n(x) 


Ra, of a 


(1) 
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像 ， 即 SDO/AIA 图 像 ，7 表示 潜在 的 清晰 图 像 ，k 表 示 
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数字 化 研究 与 应 用 示 


ASI PSF, n(x) 表示 在 每 个 像素 x=(xy) 处 相同 且 独 立 分 布 的 噪声 项 。 
根据 求解 方法 的 不 同 ， 盲 退 卷 积 算法 大 致 可 分 为 两 类 : 一 类 是 通过 引入 额外 约束 条 件 同时 对 目标 
图 像 和 PSF 估计 进行 交替 友 代 ， 并 基于 最 大 后 验 概率 估计 (maximum a posterior, MAP) 完成 育 退 卷 积 
的 算法 ， 另 一 类 是 先 估计 模糊 核 ， 再 退 卷 积 得 到 清晰 图 像 。 第 一 类 方法 对 噪声 有 一 定 的 抑制 效果 ， 但 
涉及 目标 图 像 和 PSF 估计 的 多 次 迭代 计算 ,使 得 估计 过 程 非常 耗 时 ， 且 严重 依赖 于 清晰 图 像 和 PSF 的 
先 验 特性 。 对 于 SDO/AIA 图 像 而 言 ， 很 难 获 得 清晰 图 像 的 统计 先 验 信息 。 第 二 类 方法 以 基于 显著 性 边 
缘 中 和 基于 功率 谱 的 方法 为 主 ， 基 于 边缘 的 方法 是 从 原 图 显著 边缘 估计 模糊 核 ， 这 一 方法 很 大 程度 
上 取决 于 边缘 选择 ， 更 适用 于 运动 模糊 估计 。SDO/AIA 图 像 结构 模糊 ， 无 法 较 好 的 选择 边缘 。 而 基于 
功率 谱 的 方法 直接 从 输入 图 像 功率 谱 中 恢复 出 PSF 的 功率 谱 ， 再 使 用 相位 检索 算法 求解 PSF 的 相位 。 
因此 , 本 文 研 究 了 一 种 利用 模糊 图 像 功率 谱 估 计 PSF 的 单 帧 图 像 盲 退 卷 积 算法 用 于 太阳 SDO/AIA 
图 像 增强 。 首 先 对 数据 预 处 理 (包括 去 噪 和 灰 度 值 变 换 处 理 )， 从 AIA 中 估计 出 PSF 的 功率 谱 ; 再 用 相 
位 检索 算法 恢复 出 PSF， 为 了 使 该 相位 检索 算法 对 噪声 更 鲁 棒 ， 额 外 引入 了 权重 参数 &， 限 制 使 用 PSF 
功率 谱 对 PSF 相位 估计 的 硬性 约束 。 并 对 该 算法 随机 初始 化 求解 多 次 , 避免 陷入 局 部 极 小 值 。 计 算 PSF 
所 有 可 能 解 之 间 的 L, 范 数 距离 ， 通过 多 次 随机 初始 化 计算 PSF， 如 果 可 能 解 之 间 的 LL 距离 越 小 ， 可 以 
认为 是 最 接近 真 解 的 PSF， 最 后 利用 该 PSF 于 式 (1) 退 卷 积 出 目标 图 像 。 为 进一步 减少 噪声 影响 ， 将 
竺 增强 图 像 分 隔 为 四 幅 相 邻 局 部 子 图 ， 并 按 上 述 方法 估计 每 幅 子 图 的 PSFE， 对 四 个 PSF 进行 归 一 化 和 
S 求 取 均 值得 到 恢复 整 幅 图 像 的 PSF。 实 验 分 析 使 用 轮廓 切片 曲线 、 功 率 谱 分 析 以 及 PSE 分 析 以 及 和 其 
No 他 较为 先进 的 的 增强 算法 结果 对 比分 析 ， 对 重建 出 的 目标 图 像 细节 结构 进行 客观 分 析 ， 将 高 分 辩 率 日 
© KB AL(The High-Resolution Coronal Imager, Hi-C) 图 像 功 率 谱 假设 模型 验证 和 精细 结构 增强 效果 对 
比 。 结 果 显 示 本 文 方法 对 于 SDO/AIA 193 A, 211 Å, 171Å 波段 结构 的 增强 和 复原 有 较 好 的 效果 ， 
对 噪声 鲁 棒 。 此 外 ， 由 于 该 方法 不 涉及 目标 像 的 反复 重建 ， 在 图 像 恢复 速度 上 也 有 一 定 优势 。 


Pa 


= 1SDO/AIA 数据 


N SDO/AIA"” 主要 目标 是 了 解 太 阳 能 量 存储 与 释放 机 制 ， 确 定 太阳 如 何以 及 为 何 发 生变 化 。 通 过 研 
三 | 究 太 阳 周围 等 离子 体 的 相互 作用 ， 以 增强 我 们 对 太阳 如 何 影响 地 球 大 气 和 社会 生产 的 理解 。SDO/AIA 
>: 由 四 个 望远镜 组 成 , 可 同时 提供 十 个 波段 的 全 日 面 像 , 包括 7 个 极 紫外 波段 (94A、131A、171A、193 A, 
.二 211 A, 304 A, 335 A) 和 2 个 紫外 波段 (1600 A, 1700 A) 和 1 个 可 见 光 波段 (4500 A)。 观 测 范围 包括 太 
二 阳 过 渡 区 和 日 园 上 最 高 可 达 0.5Re 的 日 晃 ， 其 空间 分 辨 率 为 1.5-arcsec， 时 间 分 辨 率 为 128， 视 场 达到 了 
41X41 arcmin, F 20104F 2 H 11 ARs. 

EDIR A ARBA CHC) HRT ARRAN SAUL, SRR TT ET PERK (50.3 ~0.4 
arcsec)， 高 时 间 分 辩 率 (25s) ARRAY Sy xX Baa AN aR AL H BN LET RAE AN tee i 7) EA H ER 
像 。 使 许多 人 了 解 到 一 些小 尺度 日 冕 特征 和 发 生 在 过 渡 区 的 活动 现象 ， 以 及 它们 之 间 可 能 存在 的 相互 
关联 也 得 到 进一步 验证 。 

如 表 1 数据 对 比 可 知 ，Hi-C 虽然 拥有 更 高 的 时 空 分 辨 率 ， 但 其 视 场 小 、 拍 摄 时 长 较 短 ， 无 法 稳定 
长 期 的 对 日 冕 进行 观测 ， 相 较 而 言 SDO/AIA 能 够 提供 大 视 场 、 长 期 全 日 面 观测 图 像 ， 然 而 其 图 像 质 
量 低 于 Hi-C， 因 此 对 SDO/AIA 图 像 增 强 处 理 显 得 尤为 重要 。 

表 1. SDO/AIA 与 Hi-C 数据 特点 对 比 表 
Table 1. Comparison of SDO / AIA and Hi-C data characteristics 


Data Properties 


Data 


Data Acquisition Field of view Image Size Spatial Resolution Tempotal Channel (A) 
Time (arcmin) (pixel) (arcsec) Resolution (s) 
SDO/AIA 2010.02.11-now 41x41 4096x4096 1.5 12 94, 131, 171, 193, 211, 
304, 335, 1600, 1700, 4500 
Hi-C 2012.07.11 6.66x7.03 <3880x4096 0.3 ~0.4 z5 193 


18:52:10-18:57:49 
本 工作 采用 SDO/AIA ŽE HE FY DA WL SS Bl) PAE BEE BS A SE 2012 4E 7 H 11 H 18:53:32 时 
刻 193A 波段 局 部 像 ， 和 观测 活动 区 日 冕 的 2014 F 2 H 2A 23:40:11 的 211A 波段 局 部 像 ， 这 两 组 数 
据 由 同一 望远镜 拍摄 。 此 外 ， 实 验 中 还 包括 不 同 望 远 镜 拍摄 的 观测 宁静 区 日 时 的 2014 年 10 月 4 日 
10:13:11 的 171A 波段 局 部 像 。 其 中 193A 波段 实验 中 ， 也 采用 了 活动 区 域 11520 的 2012 年 7 月 11 
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H 18:53:27 时 刻 Hi-C 图 像 作 为 清晰 图 像 功率 谱 假 设 验 证 和 实验 对 比 ， 除 此 之 外 ， 实 验 对 比 图 还 包括 
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文献 [3] 估 计 的 PSF 退 卷 积 图 。 数 据 预 处 理 包括 : 根据 SDO/AIA 的 数据 的 高 动态 特点 ， 对 图 像 对 数 变 
换 ， 这 一 步 能 够 有 效 降低 图 像 的 动态 范围 ， 然 后 使 用 中 值 滤 波 对 图 像 去 噪 处理， 该 算法 在 去 除 部 分 吕 
声 的 同时 不 会 平滑 过 渡 ， 能 有 效 保持 图 像 结构 。 


2 盲 退 卷 积 算法 理论 


本 文 方法 分 为 三 部 分 ， 自 然 图 像 功率 谱 假设 模型 ，Hi-C 功率 谱 分 析 和 AIA 图 像 PSF 求解 。 
3.1 自然 图 像 功率 谱 假设 模型 
Amit 等 人 研究 表明 ， 清 晰 图 像 的 功率 谱 (或 自 相 关 ) 是 各 向 同性 的 ， 模 糊 过 程 中 图 像 的 中 、 高 频 
分 量 降低 ， 对 模糊 图 像 不 同方 向 使 用 微分 滤波 后 再 计算 自 相 关 ， 当 模糊 方向 与 微分 方向 一 臻 时， 损失 
最 大 ， 能 够 认为 此 时 能 量 绝对 值 之 和 最 小 的 方向 为 模糊 方向 。 根 据 以 上 结论 ， 我 们 对 清晰 图 像 和 模糊 
图 像 分 别 进行 不 同 角度 (9:1 180") 的 传 里 时 中 心 投影 变换 ， 得 到 不 同 投影 角度 对 应 的 一 维 切片 后 ， 
再 对 这 些 切片 一 维 微分 滤波 ， 最 后 计算 其 自 相 关 即 可 得 到 : 
(Rav i (x)= cB(x) (2) 


(Rasen) (X) * CoLBty ce () (3) 
HP, d 为 一 维 微分 滤波 器 ，c 表示 一 个 随 着 角度 变化 的 乘 性 因子 ， 嘲 表示 对 图 像 自 相 关 计 算 ，P 表 
示 投 影 算 子 ，5 表示 冲 激 函数 ， 即 
sf A =(00) 


0 otherwise 


3.2 Hi-C 功率 谱 分 析 

如 图 1 所 示 ， 对 Hi-C 图 的 功率 谱 特 性 进行 分 析 。 将 图 Hi-C(a) 的 一 个 局 部 子 图 的 况 4ww,w 绘制 如 图 
(b), 图 (b) 将 45 .90 135 180° 方向 功率 谱 截 面 绘制 如 图 (c), 将 图 (c) 中 四 条 曲线 对 应 数值 相 除 得 到 结果 图 
(d)。 通 过 观察 图 (c)， 清 晰 图 像 (a) 经 过 (2) 式 计算 得 到 功率 谱 并 不 是 一 个 精确 的 8 函数 ，(d) 图 可 以 看 出 
不 同 角度 的 功率 谱 相 除 得 到 不 同 常数 ，ce 随 着 投影 角度 的 变化 而 变化 , 即 太 阳 清 晰 图 像 Hi-C 功率 谱 随 
着 6 不 同 有 不 同 乘 性 因子 的 变化 量 。 进 而 证 明 (2) (3) 式 用 于 AIA 图 像 PSF 估计 是 可 靠 的 。 图 中 多 个 尖 
峰 是 除数 存在 0 值 造成 的 。 

为 了 不 增加 网 像 的 总 能 量 ， 滤 波 器 ad 均值 为 0， 在 进行 a*Pe(B) 计算 时 ， 会 导致 Po(B) 的 均值 缺失 
(为 0) ， 但 图 像 本 身 功 率 谱 并 非 为 0， 所 以 在 计算 过 程 中 ， 引 入 参数 mm 来 代 蔡 OA. HR o (x) 的 均 
值 使 用 m/co RE E f(x) = Rar (Xx)， 由 (3) 式 可 以 得 到 : 


Rew (X)=(fo(x) + )/co (4) 


Slice of Autocorrelation('d) Stic Division 


aged 


8 
8 
B 


Or 02 
Frequency(t/pixel) 


(c) 功 率 谱 切片 图 (d) 不 同 切片 相 除 
图 1. Hi-C 功率 谱 模型 说 明 
Figure 1. Explanation of Hi-C power spectrum model 


(a)Hi-C 


3.3 PSF 估计 
(4) 式 是 一 个 欠 定 方程 ， 为 了 成 功 恢复 Rs,0) (Xx) ， 需 要 引入 额外 约束 条 件 估计 出 ce 和 mm 。 
3.3. 1 估计 Ce 


采用 关于 PSF 的 三 个 假设 "对 式 (4) 进行 等 价 计算 : 


1， 相 机 曝光 ， 光子 在 传感器 上 累积 图 像 ， 由 于 光子 能 是非 负 ， 故 PSF 是 非 负 的 ， 即 20， 所 以 关于 


的 投影 Po(k) 和 自 相关 9 都 是 非 负 的， 根据 (4) 式 能 够 得 到 态 (x)+m>0。 


2. 相机 曝光 时 移动 有 限 , 故 模糊 过 程 有 限 , 上 支持 域 紧凑 , 所 以 对 于 每 个 角度 的 投影 的 自 相 关 Re R 


Arm leg . 
Co ， 故 mm 满足 me =—Fo (So) 。 
0 otherwise 
3， 理 论 上 ， oe 
ao(x)dx=1。 根 据 式 (4 得到，cs = f(A (x) +m jax (BI ET BOI k 14) 。 
R (x)= fy(x)— fls) _ fy(x)— Als) 
J (fo (x) +m, x J (fe (x)- fa (50) x 


说 ， 都 存在 一 个 % 1849 R (x)= 


3. 3. 2 估计 支持 范围 5。 
通过 在 PSF 估计 和 8% 估计 之 间 进 行 最 大 期 望 (Expectation-maximization EM) 迭代 ， 从 而 估计 出 


和 PSF. EM 算法 的 每 一 次 迭代 都 是 对 PSF 的 再 一 次 修正 ， 能 够 更 准确 估计 PSF， 使 其 在 讶 复原 工作 


中 有 效 抑制 振 铃 和 伪 影 

(1) PAE WML sP =arg min, 万 ( 习 ， 之 后 每 次 更 新 为 中 ; 

(2) AEH st DEH RY) (x) ， 再 根据 相位 检索 算法 (算法 2) 迭代 出 当前 满足 最 大 期 望 值 的 Ko 
(3) 根据 第 (2) 步 的 ko 更 新 出 sf 用 于 第 寺 1 次 迭代 ; 

(4) (2) (3) 两 步 反 复 迭 代 直 到 满足 条 件 信 ”>0.10max (Ri ) 时 则 停止 迁 代 。 


图 2 展示 了 支持 域 壕 代 过 程 (a) 和 根据 最 终 5 过 代 结果 恢复 的 PSF 一 维 自 相关 (b)。(a) 中 蓝 色 曲线 
表示 初始 化 的 % ， 根 据 EM 迭代 更 新 % (红色 曲线 ) ， 达 到 终止 条 件 后 % (绿色 曲线 ) 的 数值 用 于 作为 (b) 


图 横 轴 的 支持 域 数值 。 


| 


| 


Angle 8 Final s, 
(a) Sp ERRE (b) PSF 功率 谱 
图 2. Sy 迭代 路 径 (a) 和 一 维 PSF 功率 谱 (b) 
Figure 2. So Iterative path (a) and one dimensional PSF power spectrum (b) 
算法 1 (求解 上 ) 的 步骤 可 归纳 为 如 下 : 
输入 : 模糊 图 像 8; 


输出 (Enk): 
Stepl. f(x)= R08) ; 
Step2. 初始 化 s =arg min, fy (x) ; 


Step3. 外 循环 开始 ， 令 次 数 t=1; 
Step4. 内 循环 开始 : 


Step4. 1. 使 用 给 定 的 %s， 根 据 式 (5) 计算 模糊 核 功率 谱 四 (“表示 对 信号 做 储 里 叶 变 换 ) ; 


Iteration process 
Angle 8 


Step4. 2. 通过 算法 2 T+ EH PRK k; 


Step4. 3. EM ERR E 5 ; 


Step5. 7# 4 So <0. ax( 


Rp o ) 时 内 循环 结束 ; 


Step6. t=T (T 表示 剪 切 图 像 张 数 ) 时 外 循环 结 
图 3 是 本 文 算法 流程 ， 首 先 对 原 图 (9) 去 品 和 对 数 变化 得 到 (pD)， 计 算 (b) 的 一 维 自 相关 (c)， 根 据 (c) 
二 维 功率 谱 (f)， 相 位 提取 算法 计算 


ARH REER so (d); ea AEH IX so BEME PSF 一 维 


出 相位 后 恢复 出 PSF(g)， 最 后 非 盲 退 卷 积 得 到 目标 图 像 (h)。 


(a) EWI EI B 


(b) log, (B) (e) Rasps) 


功率 谱 (e)， A] 原 为 


(e) Rew) 


图 3. 本 文 算法 的 执行 过 程 :输入 AIA 观测 的 


后 的 一 维 自 相关 投影 (e); 恢复 后 的 二 维 模糊 核 功 率 谱 (f) ;相位 恢复 后 的 模糊 核 (g) ;本 文 算法 增强 后 的 太阳 图 (h); 


Figure 3. The process of the algorithm in this paper: input the AIA local image (a); the logarithmic enhancement image(b); 


c [ef Co) 模糊 核 k 


(h) 本 文 算法 复原 图 7 


局 部 太阳 模糊 图 (a) ;进行 对 数 增强 后 的 图 (b); 一 维 自 相关 投影 (c); 56 迭代 路 径 (d); Sy 优化 


one-dimensional autocorrelation projection (c); Sg iterative path (d); Sg optimized one-dimensional autocorrelation projection (E); 


recovered two-dimensional PSF power spectrum (f); phase recovered PSF (g); the final enhanced image (h); 


3.4 PSF 的 相位 恢复 


RE 之 后 ， 使 用 相位 提取 算法 ， 


模糊 核 E 。 相 位 提取 算法 是 os. 该 算法 的 解 存 在 平凡 模糊 的 问题 ， 通 党 


过 添加 约束 条 件 求解 模糊 核 上 的 相位 分 量 p(o) ， 最 终 得 到 
常会 收敛 到 平凡 解 或 


局 部 极 小 值 。 典 型 的 相位 恢复 方法 以 误差 下 降 算 法 CError-reduction, ER) 为 例 ， 在 迭代 过 程 中 降低 误 


差 ， 但 其 将 不 满足 约束 的 信号 设 为 0 值 的 方式 常 


SICA, Jad 


b 现 混合 输入 输出 (Hybrid in- 


put-output, HIO) 算法 5 改进 了 此 问题 ， 对 解 设置 了 一 个 小 扰动 ， 可 使 得 迭代 的 时 候 跳出 停滞 现象 快 
速 收 你 。Lukeb5 提 出 了 一 种 叫做 松弛 平均 交替 反射 (Relaxed Averaged Alternating Reflection, RAAR) 


的 算法 ， 该 方法 能 够 使 得 收敛 更 快 更 平稳 ， 
出 的 信号 完全 满足 输入 的 功率 谱 , 当 输 入 信号 含有 噪声 的 时 候 ， 


但 该 方法 将 输入 的 功率 谱 作为 一 种 硬性 约束 ， 即 需 
会 使 得 收敛 震荡 ,和 鲁 棒 性 欠 佳 .SDO/AIA 


# 要 恢复 


图 像 噪声 复杂 ， 并 不 能 在 数据 预 处 理 阶段 完全 消除 ， 故 本 文 算 法 放宽 了 模糊 核 功率 谱 的 硬性 约束 要 求 


此 随机 过 程 ， 可 以 避免 收敛 到 局 部 极 小 值 。 其 次 ， 还 引入 了 HIO 算法 中 对 解 的 扰动 方法 〈 即 算法 2 中 


(通过 调整 权重 参数 a )， 可 提高 算法 对 噪声 的 鲁 棒 性 。 初 始 阶段 采用 随机 初始 化 相位 分 量 ， 


了 项 )， 可 以 使 得 解 跳出 停 澡 现象 快速 收敛 。 最 后 计算 所 有 可 能 解 的 L, 距 离 ， 求 出 最 优 解 。 
算法 ee ) 迭代 包括 以 下 步 又 : 


MA: p(o) =|k(a) 


输出 : k, ; 
Stepl. 外 循环 n=1; 


多 次 重复 
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Stepl. 1 随机 初始 化 相位 p(o) ([-r，r] 之 间 ); 
Stepl.2 计 算 G=F (pO; 
Stepl.3 内 循环 m=1; 
Stepl. 3.1 施加 健 里 叶 域 约束 ( 即 是 对 幅 值 P 放宽 至 Qa 倍 ): 
G,=F" ((ap +(1-a)|o)Je*** : 
Stepl. 3. 2 施加 空间 域 约 束 〈 即 PSF 非 负 性 约束 和 支持 域 约束 ): 


BG(x)+(1-28)G, (x) xeQ 
soe | 


Hot, O={x:2g, (x)a(x) <O}Ufx:x €[0.s]x[0.s]} ,为 使 m 次 迭代 中 ， 信 号 能 够 平稳 快速 的 收 
S 将 采用 HIO 算法 中 对 B 的 设置 ， 即 B=Bo+(1-Bo)(1-e™”), 本 实验 中 B=0.75, 当 出 现 xe9 


时 ， 为 G 赋予 一 个 较 小 的 值 BG(x)+(1-2B)G,(x)， 以 便 能 够 逃 出 迭代 过 程 中 的 停 灌 现象 ; 
Stepl.4 m=M u 时 ， 内 循环 结束 ; 


Stepl.5 复原 五 Ta | ， 其 中 Q={x:G,(x)< UL x: xe 0,4 x[ 0.9} i 


Step2 n=Noue H SAAR, k, 两 两 计算 L， 距 离 ， 得 到 最 优 解 ; 


算法 2 包括 两 个 循环 结构 ， 设 置 外 循环 为 了 避免 解 陷 入 局 部 极 小 值 ， 当 初始 化 值 越 接近 真 值 ， 越 


容易 收敛 到 全 局 最 小 值 。 多 次 相位 随机 初始 化 求解 E(x) ， 最 终 通 过 求 不 同 解 之 间 的 范 数 距离 ， 范 数 距 
离 最 小 的 解 可 认为 是 最 优 解 K(x) ， 设 置 内 循环 能 够 使 得 算法 平稳 快速 的 收 剑 ， 在 迭代 求解 过 程 中 交 霸 
引入 频 域 约束 ( 幅 值 约束 ) 和 空间 域 约束 (PSF 的 非 负 性 约束 和 支持 域 约束 )， 当 出 现 xs@ 时， 将 为 G 
赋予 一 个 较 小 的 值 BG(OJ+(-2 同 G:(a]， 以 便 能 够 逃 出 达 代 过 程 中 的 停 汪 现 象 。 文 献 [15] 指 出 通过 逐 


步 增 大 有 值 能 够 引导 解 稳定 收敛 。 故 本 文通 过 内 循环 将 B 值 设置 为 一 个 不 断 接近 1 的 动态 值 。 


如 式 (所 示 ， 已 知 模糊 核 上 和 模糊 图 像 B ， 求解 目标 清晰 图 像 的 方法 叫做 非 盲 退 卷 积 算法 。 实 验 


使 用 文献 [16] 中 交 禁 最 小 化 的 方法 ， 即 半 二 次 分 裂 求 解 方 法 来 求解 目标 图 像 7。 


3 实验 结果 分 析 


为 了 评估 该 方法 的 有 效 性 和 准确 性 ， 本 文 实验 中 对 AIA193A，AIA211A，AIA171A 三 个 波段 不 同 


区 域 进行 增强 ， 图 中 AIA 图 像 均 通过 对 数 显示 。 图 4 中 三 个 波段 在 增强 后 ， 细 节 结 构 更 锐 ， 尤 其 冕 环 
之 间 的 相对 模糊 结构 也 在 增强 后 较 好 的 分 离开 来 ， 边 界 也 更 为 清晰 。 如 第 三 列 本 文 算法 估计 出 的 PSF 
尺寸 大 小 可 知 ，193 和 211 波段 PSF 较 171 波段 更 大 , 这 说 明 前 两 个 波段 图 像 模糊 程度 更 大 ,由 于 193 
波段 和 211 波段 图 像 主 要 是 拍摄 太阳 活动 区 得 到 的 ， 活 动 区 更 大 的 动态 范围 和 冕 环 相对 望远镜 的 运动 
造成 模糊 程度 更 高 。 在 第 四 列 功率 谱 图 中 ， 三 个 波段 图 像 增 强 后 的 功率 谱 曲 线 相对 AIA 图 像 曲线 在 中 
频 处 都 有 一 定 程度 的 提高 。 中 频 的 提高 在 增强 图 像 中 表现 为 已 有 结构 边界 更 锐 ， 以 及 由 于 模糊 没 观 察 
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(a)AIA (b)Enhanced images (c)PSF (d)Power spectrum 
4.AIA 不 同 波段 增强 图 
Figure 4. Enhancement images of AIA in different bands 


为 了 更 好 的 对 比 细节 结构 的 增强 情况 , 我 们 在 图 5 中 引入 了 结构 的 轮廓 切片 图 , 文献 [6] 根 据 ATA 
望远镜 的 理论 估计 出 不 同 波段 PSF 用 于 增强 AIA 图 像 ， 太 阳 一 个 波段 的 所 有 图 像 均 使 用 同一 个 PSF 
退 卷 积 得 到 ， 这 是 不 准确 的 ， 图 像 模 糊 原 因 也 只 考虑 了 望远镜 设备 的 因素 。 所 以 恢复 结果 不 佳 ( 如 图 
5 第 二 列 )， 图 5 中 ， 我 们 对 不 同 波段 不 同 的 区 域 分 别 增强 ， 并 将 本 文 估计 的 PSF 与 文献 [6]PSF 作对 
比 。 如 193 波段 的 两 张 局 部 图 ， 本 文 估计 的 PSF 尺寸 相当 ,但 形状 不 同 ， 对 于 结构 A 和 结构 BB, RR 
被 较 好 的 分 离开 来 ， 如 轮廓 切片 图 所 示 。 同 样 ， 在 211 波段 和 171 波段 增强 图 中 ，C、D、E 结构 的 模 
糊 方向 与 PSF 形状 一 致 ， 不 同 的 区 域 ， 太 阳 的 活动 不 一 致 ， 冕 环 位置 和 形状 的 变化 ， 相 互 干扰 是 不 可 
控 的 。 这 也 是 必须 分 区 域 单独 估计 PSF 退 卷 积 图 像 的 原因 。 


Normalized Intensity 
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(a)AIA (b) Method of reference [6] (c)Our method (d) Contour slice map 


图 5. 不 同 波段 不 同 区 域 细节 增强 图 


Figure 5. Detail enhancement images of different bands and regions 
如 图 6 所 示 ， 通 过 引入 更 高 分 辨 率 的 Hi-C 图 像 进行 细节 结构 对 比 ， 评 估 本 文 增强 方法 的 准确 性 。 
为 了 便于 与 Hi-C 图 对 比 ， 图 (a)(b) 经 过 放大 和 对 数 处 理 ， 观 察 三 者 的 能 量 图 可 知 ，AIA 图 像 细节 区 域 
颜色 相近 ， 模 糊 明 显 ， 几 乎 看 不 到 细节 结构 ， 增 强 后 的 图 像 与 Hi-C 相近 。 通 过 三 者 的 功率 谱 图 (d) 可 
以 看 出 ，AIA 几乎 没有 高 频 细节 ， 中 频 部 分 也 很 弱 。 经 过 本 文 方法 增强 后 的 图 像 ， 中 高 频 都 有 一 定 的 


提高 ， 略 弱 于 Hi-C 图 像 ， 说 明 恢 复 细节 结构 的 数量 是 少 于 Hi-C 


图 像 的 ， 中 频 处 有 一 点 略 高 于 Hi-C 


像 , 在 细节 图 (e) 中 可 观察 到 ， 本 文 方法 不 仅 能 较 好 的 分 离 园 环 的 相对 模糊 结构 ， 且 对 于 已 存在 的 金 环 ， 
能 更 好 的 保持 边缘 。 细 节 结 构 边 缘 较 Hi-C 锐 ， 所 以 中 频 曲 线 略 高 于 Hi-C。 图 中 的 A、B、C 结构 都 是 


由 于 模糊 无 法 分 辨 的 结构 ， 增 强 后 ,与 Hi-C 对 比 可 知 ， 结 构 保 持 稳定 ， 且 恢复 情况 较为 准确 。 如 对 应 


的 轮廓 切片 图 (DA 结构 所 示 ，AIA 中 表现 为 一 个 锡 环 的 结构 ， 成 功 的 分 离 为 两 条 金 环 ,与 Hi-C 保持 一 


致 。 
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(£) Contour slice map 
图 6. 与 Hi-C 细节 对 比 图 
Figure 6. Detail comparison with Hi-C 


图 7 给 出 了 经 典 的 基于 直方 图 的 BPDHE 算法 ， 以 及 目前 较为 先进 的 EFF, MF, rR 算法 的 实验 结 
果 ， 第 一 行 局 部 放大 细节 中 能 看 出 ，BPDHE 方法 使 得 图 像 局 部 过 度 增 强 且 细节 模糊 ， 本 文 算法 的 增强 
图 对 细节 增强 效果 更 好 ， 边 缘 结 构 也 较为 稳定 。 第 二 行为 各 个 算法 对 应 的 直方 图 ， 根 据 直 方 图 像素 值 
的 分 布 情况 可 看 出 ，ME 算法 的 可 视 化 效果 最 好 。 对 比 其 他 增强 算法 ， 本 文 算法 在 保证 对 比 度 和 亮度 
增强 的 同时 ,能 复原 出 更 精细 的 细节 结构 。 如 表 二 所 示 ， BRA DT AS TR IE ad oh ET be A I o 
这 两 者 可 有 效 的 反映 图 像 的 对 比 度 和 亮度 信息 。 其 值 越 大 , 表示 图 像 信 息 丰 富 。 亮度 误 差 [] (Lightness 
Order Error, LOE) 能 够 客观 的 衡量 图 像 的 亮度 失真 情况 ， 值 越 小 代表 亮度 越 好 。SSIM (structural 
similarity) 可 用 来 衡量 图 像 恢复 质量 ， 是 图 像 亮 度 ， 对 比 度 和 结构 对 比 三 者 均 评 估 的 指标 。 其 值 越 接 
近 于 1， 图 像 质量 越 好 。 表 中 的 方法 均 与 Hi-C 图 像 作对 比 计算 SSIM。 相 较 而 言 ， 虽 然 本 文 算 法 运行 
时 间 较 长 ， 但 对 细节 结构 的 增强 更 为 准确 。 


$ $- 


到 7. 不 同 增强 算法 对 比 


Figure7.Comparison of different enhancement algorithms 


表 二 . 图 像 增强 测量 参数 


Table2.Image enhancement measurement parameters 


AIA BPDHE EFF MF rR Ours 
Standard deviation 647.7 953. 1 1768 1974. 6 2270. 1 2416.9 
Entropy 2.679 6. 428 7. 223 7. 653 7. 306 11. 873 
LOE 871 926 374 326 344 349 
SSIM 0. 2708 0. 210 0. 303 0. 435 0. 326 0. 716 
Running time(s) (400x400) 一 1. 033 1. 134 1. 332 21. 151 13. 718 


44 论 


本 文 提出 了 一 种 基于 盲 退 卷 积 的 太阳 AIA 图 像 质量 增强 方法 , 基于 图 像 功 率 谱 特性 的 PSF 估计 模 
型 ， 从 ATA 图 中 估计 出 PSF， 最 后 退 卷 积 出 增强 图 像 。 文 中 采用 轮廓 切片 分 析 评 佑 方法 准确 性 ， 功 率 
谱 分 析 评 佑 方法 有 效 性 ， 不 同 区 域 PSF 分 析 评 估 方 法 的 合理 性 。 实 验 结果 表明 ， 本 文 方 法 对 
SDO/AIA193A,171A,211A 波段 图 像 增强 效果 较 好 ， 尤 其 纹理 细节 ， 边 缘 结构 的 增强 效果 最 好 。 其 次 ， 
本 文 算法 也 有 一 定 的 噪声 抑制 效果 ， 得 益 于 PSF 相位 估计 过 程 的 松弛 设计 ， 恢 复 新 细节 结构 也 较为 准 
确 ， 因 为 这 三 个 波段 的 图 像 噪声 较 小 ， 日 冕 结构 特征 也 更 丰富 和 明显 ， 能 够 更 为 准确 的 估计 出 模糊 核 。 
对 比 其 他 图 像 增 强 方法 ， 本 文 方法 的 实时 性 较 差 ，PSF 估计 和 退 卷 积 的 最 优化 迭代 过 程 是 主要 的 耗 时 
原因 。 将 来 可 以 通过 提高 计算 设备 的 性 能 和 采用 较为 快速 且 有 效 的 最 优化 方法 来 提高 算法 的 实时 性 。 
在 ATA 的 其 他 波段 图 像 中 ， 同 样 可 以 观测 日 冕 的 SDO/AIA94A, 131A, 335A 波段 图 像 噪声 很 大 ， 异 党 
点 较 多 ， 结 构 相 对 模糊 导致 增强 效果 不 佳 ， 将 来 能 够 通过 研究 更 好 的 去 噪 算法 和 对 噪声 更 加 和 鲁 棒 的 退 


卷 积 算 法 对 其 他 波段 有 效 增 强 。 


oo 
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AJA Images Enhancement Method Based on Blind Deconvolution 
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Abstract: There are various physical activities in the solar image, such as different scales, brightness 
and structure. Due to the high dynamic activity of the solar corona and sensor equipment, the image qual- 
ity of the solar image is poor. According to the characteristics of large dynamic range, large noise and rel- 
atively fuzzy structure of different waveband data structures captured by the atmospheric imaging asset 
(AIA) of the solar dynamic Observatory (SDO), an effective SDO / AIA image enhancement method based 
on blind deconvolution is proposed in this paper. Firstly, the SDO / AIA image is denoised and the dy- 
namic range is reduced. Based on the distribution assumption of image power spectrum, the point spread 
function (PSF) is estimated from the original image Then, phase retrieval algorithm is used to recover the 
phase of PSF, and then deconvolute the high-quality target image; finally, the enhancement results are an- 
alyzed quantitatively and qualitatively through slice analysis, power spectrum analysis and PSF analysis. 
The experimental results show that, compared with the existing image enhancement methods, the pro- 
posed method can effectively enhance the detail structure of the solar corona image, and restore the struc- 


ture that cannot be identified due to the blur of the original image. 
Key words: solar image; SDO/AIA; power spectrum; blind deconvolution; phase retrieval 
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